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ABSTRACT 



We use Hubble Space Telescope Fine Guidance Sensor astrometry and high- 
cadence radial velocities for HD 136118 from the HET with archival data from 
Lick to determine the complete set of orbital parameters for HD 136118 b. We 
find an orbital inclination for the candidate exoplanet of % = 163.1° ± 3.0°. This 
establishes the actual mass of the object, M b = 42+JgMj, in contrast to the 
minimum mass determined from the radial velocity data only, M b sini ~ 12 Mj. 
Therefore, the low-mass companion to HD 136118 is now identified as a likely 
brown dwarf residing in the "brown dwarf desert" . 
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1. INTRODUCTION 

Among the hundreds of exoplanets detected so far, only a few of them have their actual 
mass known. The widely used Doppler spectroscopy technique yields the radial component 
of the stellar perturbation velocity only. Consequently, the inclination of the orbital plane 
is unknown and only the minimum mass of a companion may be determined. To obtain a 
companion's true mass it is necessary to make use of additional techniques. 

The first precise determination of an exoplanet mass was made for a transiting system 



(Henry et ai. 2000). However, transits are observed to occur only for systems that are 
oriented edge-on, and they have only a reasonable probability of occurrence for close-in 
planets (semimajor axes less than about 0.1AU). Another way to determine the orbital 
inclination of an unseen companion is by measuring the stellar reflex motion astrometrically. 



The first astrometrically determined mass of an exoplanet by Benedict et ai. (2002b) was 
possible thanks to the high precision of the Fine Guidance Sensor (FGS) instrument on the 
Hubble Space Telescope (HST). The FGS provides per observation precisions of better than 
lmas for small angle relative astrometry. This unique capability enables the detection of 
stellar perturbations due to planetary mass companions in wide orbits. 

We were granted observing time with the HST to measure the perturbation and deter- 
mine the true mass of HD 136118 b, which is an exoplanet candidate found by radial velocity 



measurements (Fischer et ai. 2002). To supplement previously published data and pro- 
vide better constrain on the companion's spectroscopic orbital parameters we also obtained 
high-cadence radial velocity measurements with the Hobby- Eberly Telescope (HET). In this 
paper we present the results of our analysis and it is arranged as follows, in Section [2] we 
review stellar properties for HD 136118. In Section[3]we discuss the instrumental set up and 
data reduction for both spectroscopy and astrometry. In Section [4] we describe the orbital 
model used to analyze this data. In Section [5] we discuss the strategy employed to obtain 
the system parameters and present our results. In Section [6] we summarize and discuss the 
consequences of our results. 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Tele- 
scope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. Based on observations obtained with the Hobby-Eberly Telescope, 
which is a joint project of the University of Texas at Austin, the Pennsylvania State University, Stanford 
University, Ludwig-Maximilians-Universitt Mnchen, and Georg-August-Universitat Gottingcn. 
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STELLAR PROPERTIES 



HD 136118 (=HIP 74948) is a V = 6.93, F9V star with roughly solar photospheric 
abundances (Gonzalez & Laws 1 |2007| ). Table g] summarizes its observed properties given in 



the literature. Observations of Ca II H and K lines (Fischer et ai. 2002) indicate modest 



chromospheric activity for this star, therefore not many spots should be expected. Using 



the Saar & Donahue (1997) relationship for the spot radial velocity amplitude versus filling 



factor: As = 6.5 fg 1 v sin i, where fs is spot filling factor in percent, vsmi is the projected 
velocity in km/s and As is the spot radial velocity amplitude in m/s. For the radial velocity 



amplitude of 211 m/s, and the measured velocity wsini = 7 ±0.5 km/s (Butler et ai. 2006), 
we obtain a spot filling factor of about 6%, i.e. about 60 milimag variations. As shown in 



Benedict et ai. (1998) and Nelan et ai. (2007) the Fine Guidance Sensor itself is a millimag 



precision photometer. The variations we see in HD 136118 over 700 days are on order 4 
parts per 1000, about 4 millimags, as shown on Fig [TJ This implies only small variations 
in spectral line shapes, which typically introduces noise on the order of 5 — 10 m/s in the 
velocities. 



Fischer et ai. (2002) have also provided evidence that the stellar rotation period is 



about 12.2 d. Given a stellar radius of R = 1.58 ± 0.11 R & (Allende Prieto & Lambert 
1999) we calculate the maximum rotation speed at the stellar equator v max = 6.5 ±0.2 km/s. 
The measured projected velocity v sinz = 7 ± 0.5 km/s is then consistent with the maximum 
speed. This suggests a very high inclination of the spin axis. Therefore, if the whole system 
follows the same angular momentum orientation as that of the star, then the companion's 
orbit would be close to an edge-on orientation with respect to our line of sight. 



OBSERVATIONS AND DATA REDUCTION 



3.1. HET Spectroscopy Data 



Spectroscopic observations were carried out with the High Resolution Spectrograph 



(HRS; Tull 1998) on the HET at the McDonald Observatory using the iodine absorption 
cell method (Butler et ai. 1996). Our observations include a total of 168 high resolution 



spectra which were obtained between UT dates 2005 December 4 and 2008 May 20. Multiple 
observations were taken most nights, so the velocities obtained on the same night were 
combined, producing individual measurements of the stellar Radial Velocity (RV) at 61 
different epochs. 



The spectrograph was used in the R = 60, 000 mode with a 316 line/mm echelle grating. 
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The position of the cross dispersion grating was chosen so that the central wavelength of 
the order that fell in the break between the two CCD chips was 5936 A. A temperature 
controlled cell containing low pressure iodine (I2) gas was placed in front of the spectrograph 
slit entrance during all the exposures. The absorption of light by the I2 gas produces a 
set of well known spectral features imprinted at the same time as the stellar spectrum. 
However they are produced with no wavelength shift with respect to the observatory frame. 
This provides a much better reference for wavelength calibration and also permit us to 
characterize the instrumental profile with great accuracy. 

The exposure times were nominally 120 s, but were increased on a few nights due to 
bad seeing conditions. In addition to the program spectra we have also obtained HD 136118 
template spectra. For these we removed the I 2 cell, the resolution was set to R = 120, 000, 
and the exposure times were 230 s. 

A detailed description of our reduction and radial velocity analysis of HET HRS data is 



given in Bean et al. (2007). This provides us with measurements of the stellar radial velocity 
relative to an arbitrary zero point obtained from individual spectra. The offset of the entire 
data set is determined simultaneously with the orbital fit, then after compensating for this 
offset, we combined our RV data with previously published velocities from Lick Observatory 
( Fischer et al. [2002 ) to produce a total data set that spans 10.3 yr. Table [2] contains reduced 



HET data for the observed epochs. 



3.2. HST Astrometry Data 



Astrometric observations were obtained with the Fine Guidance Sensor lr (FGS-lr), a 
two-axis, white-light interferometer aboard HST. Our data were obtained with the FGS in 
fringe-tracking position mode. A detailed description of this instrument is found in |Nelan 



et al. (2007). All observations were secured under 2-gyro guiding, an operational mode 



dictated by gyro failures on HST. This mode results in major constraints on HST roll angle 
and observation dates. The dates of observation, the number of measurements for each date 
(N bs) and the FGS-lr orientation angles are listed on Table [3j Our data sets span 1.8 yr, 
covering about 55% of the companion's orbital period. 

In order to obtain high astrometric precision (~ lmas) with FGS we perform the fol- 
lowing reduction procedures to remove the main sources of systematics. We apply an optical 



field angle distortion (OFAD) calibration (Whipple et al. 1995 McArthur et al. 1997); 



apply time-dependent corrections to the OFAD (McArthur et al. 1997) with an additional 



as yet unpublished improvement due to additional up-to-date M35 data that has provided 
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better models to recalibrate the telescope; and we correct for drift during each observation 



set (intra-orbit drift). These procedures are described in detail in Benedict et al. (1999 



2002c) and McArthur et al. (2001). 



The FGS provides the measurement of each star position in a serial fashion. Each date 
listed in Table|3]contains multiple measurements, alternating between the target (HD 136118) 
and the reference stars (REF-14, REF-16 and REF-17), comprising a total of about 4 mea- 
surements per star. This provides x'(t) and y'(t) coordinates at a time t of an epoch, which 
are in the HST reference frame. 



A neutral density filter (F5ND) was applied when observing HD 136118 due to its bright- 
ness. For the reference stars we used the F583W filter. 



3.3. Classification Spectra and Photometric Data 

We obtained classification spectra data for the astrometric reference stars. We include 
this in the models as bayesian a priori data to improve the accuracy of our determination 
of the orbital parameters for HD 136118 b. The spectra were obtained at the KPNO 4m 
Telescope and the photometric data at the NMSU lm telescope in May of 2006. Table [4] 
summarizes spectral and photometric information for the reference stars. 



3.4. Astrometric Model 



Given the positions (x', y') measured by FGS-lr we build a model that accounts for 
positional changes occurring systematically in all reference stars from date-to-date. This is 
accomplished by solving an overlapping plate model which includes scaling-rotation (Ci, C2, 
C3, C4) and offset (Di,D 2 ) constants, which are constrained to an arbitrary frame adopted 
as the reference (the constrained plate). The astrometric model also accounts for the time 
dependent movements of each star, given by the absolute parallax n abs and the proper motion 
components; \i a and (is, where the indices specify the direction in equatorial coordinates. 
Therefore the model is given by the standard coordinates £ and rj: 



£ = C X X' + C 2 y' + Di - PalTabs - (l a At (1) 

7] = C 3 x' + C 4 y' + D 2 - P 5 7r abs - fi s At (2) 
where P a and Ps are parallax factors obtained from a JPL Earth orbit predictor (Stan- 
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dish 1990). These are called the equations of condition, which comprises two equations for 
each star and for each epoch, giving a total of 144 equations to be solved simultaneously. 
In order to find a global solution we used a program written in the GAUSSFIT language 
( Jefferys, Fitzpatrick fc McArthur~||1988 ). Table [5] contains the resulting astrometric catalog. 



Van Leeuwen (2007) has recently presented a new reduction of Hipparcos data. This 



resulted in an absolute parallax of HD 136118 7r 6 s = 23.48 ± 0.54 mas, which differs signifi- 
cantly from 7r a f, s = 19.13 ± 0.85 mas, the prior Hipparcos determination, and from our HST 
result (7T a bs = 19.12 ±0.22 mas). When we constrained the parallax in our models to the new 
higher Hipparcos value, we found that the \ 2 °f our solution was increased by 20%, which 
offers evidence that in this case, the older Hipparcos determination was more accurate, if 
not more precise. 



4. ORBITAL MODEL 

4.1. Radial Velocity Model 

The velocity we are modeling is the radial component of the stellar orbital movement 
around the barycenter of the system, which is given by the projection of a Keplerian orbital 
velocity to observer's line of sight plus a constant velocity T. This constant in practice is 
not the actual velocity of the whole system but a parameter that absorbs all unaccounted 
for constants. Therefore we have the following equation: 



v = T + if [cos (/ + uj) + e cos u] (3) 

where to is the argument of periastron, e is the eccentricity, K is the velocity semi- 
amplitude and / is the true anomaly. The latter gives the time dependence, which is obtained 
implicitly by solving the Kepler equation: 

2 *{t-T ) = E-esmE (4) 

where To is the epoch of periastron passage, P is the orbital period and E is the eccentric 
anomaly, which is related to / by the following equation: 



/ fl + e E 

tan — = \ tan — 

2 V 1 - e 2 



(5) 
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The velocity semi-amplitude can also be written in terms of orbital elements: 

P y 1 

where i is the orbital inclination and a s is the semi-major axis of the orbit of the star. 
Using the proportionality between the masses and semi-major axes, a^M^ = a s M s , and 
Kepler's third law, we can rewrite Eq[6j 



Mi, sin i 
(M b + M s )2/3 




Note that we have introduced the indices s and b to distinguish between stellar and 
companion's parameters. Eq[7]provides a way to calculate the project minimum mass Mf, sin % 
of the companion with the assumption that there is a measurement of the stellar mass by 
other means (e.g. stellar atmospheric models). It is to be noted that this is a lower limit to 
the mass with the uncertainties mostly dominated by the determination of the stellar mass. 



4.2. Apparent Orbit Model 

Astrometric data provides positions of the parent star on the plane of the sky at different 
epochs. The high precision of FGS allows us to measure the apparent orbital movement of 
the star due to the presence of its companion. Therefore the astrometric observables to be 
modeled are the coordinates of the star apparent orbit. 

First we write the elliptical rectangular coordinates x,y, in the unit orbit, given by: 



x = (cosE — e) (8) 
y = ^/l^sinE (9) 

where e is the eccentricity and E is the eccentric anomaly. Remember that E carries the 
dependence on time through the Kepler's equation. The projection of this true orbit onto 
the plane tangent to the sky gives the coordinates Ax, Ay. This projection can be expressed 
mathematically by: 



Ax = Bx + Gy 



(10) 
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Ay = Ax + Fy (11) 
where B, A, G, F are the Thiele-Innes constants, given by: 

B = a(cosu; sinf2 + sin u; cos f2 cos i) (12) 

A = a(cosu; cosf2 — sin a; sin f2 cosi) (13) 

G = a(— sina;sinfi + cos u; cos cos 2) (14) 

F = a(— sinu; cos Q — cos to sinf2 cosi) (15) 



where a is the semi-major axis, Vt is the longitude of the ascending node, % is the 
inclination of the orbit plane to the plane tangent to the sky and to is the argument of 
periastron. It should be noted that this can be the orbit coordinates for either the parent 
star or its companion around the barycenter, depending on which semi-major axis is taken; 
the uj in the respective orbits differ by 180°. We measure the star orbit, so the coordinates 
of interest are Ax s ,Ay s , obtained by taking a = a s , which is expressed in mili-seconds of 
arc. 



4.3. Combining RV and Astrometry 

The orbital elements (P, e, lu, K) in the RV model are the same as the ones used in the 
astrometric model. However, RV data has a time baseline much longer and a sampling much 
more abundant than that of the astrometric data. Therefore, in order to obtain the unknown 
parameters a s , i and Q, we will use the RV parameters to enforce a "constraint relationship" 
between the astrometric and RV data sets. A way to constrain the orbit without including 
the time dependence is by making use of Eq|6j Rearranging the terms we have the following 



equation (Pourbaix & Jorrisen 2000): 



a s smt = (lb J 

where on the left side we have the projection of the semi- major axis of the perturbation 
orbit. The quantity a s is an astrometric observable obtained in angular measure. It may 
be converted to linear measure (AU) by the relation a[AU] = a/7r abs . The right hand side 
carries all quantities obtained from RV. 

If parameters obtained from RV analysis are assumed as constants, there will be three 



independent measurements, Eqs 10, uA\ and [161 to determine only three parameters, a s , i 
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and £7. It means that even if one has a poor astrometric coverage of the orbit, it would 
still converge to a unique solution, although uncertainties in the resulting orbital parameters 
could be reduced with additional coverage. 

The final solution is obtained by performing a simultaneous fit of Eqs[3] [v (P, T, e, u>, K, T; t)] , 



10 [Ax s (P,T,a s ,i,u,tt;t)] and 11 [Ay s (P,T,a s ,i,oj,tt;t)\, for RV and astrometry data. 



Some parameters are constrained via Eq 16 and all parameters are in some way constrained 
by their physical meaning. 



5. RESULTS 



All solutions presented below were obtained using GAUSSFIT (Jefferys, Fitzpatrick & 
McArthur 1988) programs that minimize x 2 f° r t ne equations shown in the models. As in 



Section 4.1 above index b stands for the companion and s for the parent star. 



5.1. Radial Velocity Solution 

The best fit model solution yields to orbital parameters shown on Table [6j Fig [2] shows 
HET RV data plotted together with previous published Lick data ( Fischer et al. |[2002 ). The 



fitted model is also plotted. The residuals are shown in the bottom panel. Table [7] shows the 
statistical quantities obtained from the residuals, where we used the following terminology. 
Central values: mean and median. Standard Deviation: o. Sum of squared normalized 
residuals: x 2 ■ Degrees of Freedom: DOF. Reduced chi-square: xl — xV^OF. We estimate 
these quantities for each data set separately and also for all data sets combined together. 
These allow us to evaluate the legitimacy of the fit. 



From Table [7] we note that xl f° r Lick data is close to unity, indicating satisfactory 
agreement between the dispersion of residuals and the individual errors. However this agree- 
ment between dispersion and errors is not as definitive for the HET data. We firstly check 
whether the residuals follow a Gaussian distribution and inspect the errors involved. Fig [3] 
shows the histogram of distribution of residuals for the two individual datasets separately, 
Lick and HET, and also for both datasets combined (hereafter "ALL"). We can see the 
different dispersion for each dataset. We note that each individual dataset, either HET and 
Lick, are not exactly following a normal distribution. We call attention to the fact that the 
dispersion on HET residuals is about three times larger than the error (~ 3m/s) estimated 



from previous work [e.g. Bean et al. (2007)]. This discrepancy may be identified with an 



unaccounted for systematic effect. Below we investigate the detection limits for any further 



periodic signal that could still be present in our data. 



5.2. Limits on Additional Periodic Signals in the RV Data 

The customary method for searching periodic signals in unevenly spaced data is by 



means of the Lomb-Scargle Periodogram (LSP) (Scargle 1982). Figs [4] and [5] show the LSP 
of residual RV data for two different datasets respectively: ALL and HET. The sets are 
analyzed separately because they have different errors (see Table [HI). The power in the LSP 
is weighted by the overall variance, therefore if one mixes two sets with different variances it 
would result in an overestimated power for higher levels of noise. The downside of analyzing 
datasets separately is that sampling becomes different as you have different time coverage 
and it may affect the detectability for some frequencies. 



The LSP for Lick data does not seem to show any expressive power (Fischer et al. 
2002). The LSP for HET dataset (Fig [5]) shows some peaks at the limit where False Alarm 
Probability (FAP) is as low as about 1%. The combined dataset also shows some power 
below the level of 1% FAP. This indicates that either there is still an unaccounted periodic 
signal or the sampling for those frequencies is poor. The latter may be analyzed by the 
method we describe below. 

We introduce a quantity to evaluate how much we can trust some high power found at 
a given frequency based on sampling for that frequency. We call this quantity the Amount 
of Information in the Phase Diagram (AIPD). It is defined by the following expression: 



N b 

/=-^^ln Pi /lniV 6 (17) 

i=l 

where Nj, is the number of bins in the phase diagram, pi is the probability of finding a 
data point within a given bin i, and may be calculated by Pi = rii/N, where rii is the number 
of points inside the bin i and N is the total number of data points. Note that the dependence 
on the period arises from the construction of the phase diagram. This quantity is normalized 
and therefore it varies from to 1. When I = it means that all data points are found within 
a single bin in the phase diagram and hence sampling for that frequency is very poor. From 
the other hand when I = 1 it means that probability pi is the same for every bin and data 
is equally distributed along all bins. This gives you an ideal coverage of the phase diagram. 
We call attention to the fact that AIPD does not measure the statistical significance of the 
number of data points but only the significance of how well distributed are these points in the 
phase diagram. Therefore an issue of concern is the choice of an adequate Nj,. We suggest the 
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choice is made in the same fashion as when you build a traditional histogram for inspecting 
probability distributions. Our choice of iVj, is that of N/N^ > 30 data-points-per-bin. If N 
is small such that makes Nj, < 10 then we assume Nb = 10. 



Figs [4] and [5] also show a plot of AIPD (in the plot it is multiplied by 20 and shifted 
for the sake of visualization). We note that some of the high power periods in the LSP also 
presents a decreasing on the AIPD, which means a deficit of sampling for those periods. 
This is evident for the lyr period where there is always lack of data for some part of the 
phase diagram. If you fold the lyr phase diagram twice there will still be some lack of data 
coming from the 1 yr sampling problem. This can be seen from the smaller decreases at half 
year period. Although it is smaller it may still affect the LSP. If one disregards the powers 
at periods which are close integer fractions of a year there will be no significant power left 
on the LSPs. 

However, the existence of an additional signal cannot be ruled out by looking only at 
the LSP for the following reason. If the signal comes from an orbit which follows the RV 
model (Eq[3]), then it may be considered solutions for eccentric orbits instead of strict sine 
and cosines as in the LSP. An alternative periodogram using the orbit solution is explored 
in Gregory (2007). We propose a strategy to find a possible hidden signal in our particular 



case, although it could be expanded and applied for any other system. 

Our strategy consists in making multiple attempts to model the RV, using a 2-companion 
model (linear superposition of two Keplerian orbits), and keeping the trial periods as con- 
stants in the fitting process. This process forces the minimization algorithm to search for 
the best solutions for each chosen period. This approach could result in a lower x 2 when 
including a hidden periodic component. Fig [6] shows a x 2 map over a range of periods for a 
2-companion model fitting RV HET and Lick data simultaneously. The grid point resolution 
is 1100 x 90, which means a step on the trial periods of about ~ 0.3 day for both components. 
From Fig [6] we can see a region around 255 days where we found an island of lower x 2 (darker 
regions), which indicates the presence of an additional signal. We note that the relatively 
high power peak at ~ 95 days in the LSP (see Fig [5]) is now ruled out, because any attempt 
of fitting a secondary orbit with this period results in larger x 2 ■ This approach does not 
prove the existence of another companion in the system. Rather it shows an effective way of 
finding solutions that are considerably improved by adding a weak periodic signal that could 
not be detected in the LSP. Moreover, our model uses Keplerian orbits, which look not only 
for a periodic signal but for a signal with the shape of an orbit. 

Then, we performed a refined fit using the Levemberg-Marquardt method and then a 
robust fit method for a 2-companion model with an additional signal with period of about 255 
days. This solution presents a notable improvement on the minimization of x 2 if compared 
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to the 1-companion model (see Tables [7] and [8]) . The parameters for a 2-companion model 
are shown on Table [9] From Table [8] we note an improvement on all statistical quantities. 
The a dispersion for HET data is now in agreement with that we expected. Figs [7] and [8] 
show the phase diagram of RV data and the respective component orbit model. 

Such solution suggests the presence of an additional lower mass companion to the sys- 
tem. However this would be a very eccentric orbit and considerably large Msinz planet, 
which makes us believe that such body could hardly coexist with a brown dwarf orbiting the 
same system at the distances they seem to be. In fact we have performed a stability anal- 



ysis using the MERCURY package (Chambers 1999). We have input the two companions 



around HD 136118, and considered the minimum mass, which minimizes interactions. We 
also explored the full range of inclinations for the second component. The system always 
becomes unstable for very short time scales. Besides the stability constraint, this detection 
is at the limit of our instrument, and also if we look at Fig [8] we note that only a few data 
points are contributing to form the orbit we have obtained. For these reasons we prefer 
to be cautious and take this as a "nuisance orbit" to fix an unknown source of systematic 
error present in the HET data, although the possibility of a second companion is not out 
of the question. Further re-reductions and observations will be done soon to investigate the 
origin of this signal. The parameters adopted in the following sections are those from a 
2-companion model. The modelling of the second component in the RV has a very marginal 
effect on the parameters of the astrometric detection, so it is not "polluting" our result. 



5.3. Simultaneous RV and Astrometry Solution 

The system parameters obtained in Section [5] from RV analysis are intially adopted as 
constants as we search for orbital solutions in the astrometric data. Once close to the solution 
we free all parameters and obtain the best fit model for astrometry and radial velocity 
simultaneously. We obtain a semi-major axis of the perturbation orbit a s = 1.45 ± 0.25 mas, 
an inclination % = 163.1° ± 3.0° and a longitude of the ascending node Q = 285° ± 10°. 
Fig [9] shows the reduced star positions Ax s and Ay s , versus time for HD 136118. Although 
our solutions were obtained considering each data point individually, in Fig [9j in order to 
provide the reader a better visualization to show how the fit works, we also plot normal 
points which are the median and respective standard deviation of the mean of each clump of 



data, representing 3 different epochs. These collapsed points are also shown on Fig 10 , where 



we plotted Ay s versus Ax s and the apparent orbit fit model. Fig [TT] shows the distribution 
of astrometric residuals of Ax (left panel) and Ay (right panel) for all reference stars and 
all data sets. A Gaussian fit model is superposed for comparison. The fit distribution for 
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both Ax and Ay present a maximum consistent with zero, and a FWHM of 0.87 mas and 
1.02 mas respectively. 

By determining the inclination we are able to remove the previous degeneracy on the 
mass of the companion. We calculate the actual mass by iterating Eq [7j which yields 
M b = 42±\lMj, firmly establishing HD 136118 b bloom in the brown dwarf desert. 
We also obtain the physical semi-major axis of the companion orbit, a b = 2.36 ± 0.05 AU. 
A summary containing all HD 136118 b parameters derived from the simultaneous RV and 



astrometry solution is shown on Table 10 



6. SUMMARY AND DISCUSSION 



HD 136118 is a solar type star with a brown-dwarf companion. Table 10 shows a sum- 
mary of all observed orbital elements of the system. 

We found that HD 136118 b has an orbital inclination of i — 163.1° ± 3.0°, nearly 
perpendicular to the inferred inclination of the stellar spin axis. This misalignment is an 
intriguing result since conservation of angular momentum would favor alignment between 
stellar spin and companion orbital axes, assuming both were formed in the same primordial 
cloud. 

HD 136118 b is likely a brown dwarf companion orbiting at 2.36 AU that falls in the 



driest region of the so called 'brown dwarf desert' (Grether & Lineweaver 2006). They 



showed that the frequency of companions in the stellar mass range follows a slope with 
gradient — 9.1±2.9, while in the planetary mass region, the gradient is 24.1 ±4.7. These two 
separate linear fits intersect below the abscissa at M — 43^3 Mj. Surprisingly, HD 136118 b 
mass is M b = 43^Jg Mj. Reffert & Quirrenbach (2006) measured astrometric masses for two 
exoplanet candidates HD 38529 c (M = 37l?| M Jtip ) and HD 168443 c (M = 34 ± 12M Jup ). 
Both are likely brown dwarf companions around solar type stars like HD 136118 b. These 
objects are important cases for studying the mass function at the brown dwarf mass range. 



According to the evolutionary dusty model for brown dwarfs of Baraffe et al. (2001) 
and Chabrier et al. (2000), assuming M b = 0.0411q;q}° M Q and the age of the brown dwarf 
as 5 Gyr, HD 136118 b has a temperature of about Tj, = 900 K and a radius of R b = 0.086 R & . 
If one considers the uncertainty in the age of the system (see Table [T]), this brown dwarf may 
be much younger, therefore considering the age as 1 Gyr, HD 13118 b has a temperature of 
about T b = 1200 K and a radius of R b = 0.1 R Q . These characteristics classifies HD 136118 b 
as a T dwarf. Using these values we calculate the emission and reflection spectra, and the 



flux ratio between the brown dwarf and the parent star as show in Fig 12 The flux ratio 
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increases toward the far infrared (L, M and N bands), where it can get as high as 10 



The astrometric determination of the mass of a low mass companion can decisively 
characterize it as a planet. A good illustration of this fact can be seen from the results of 
our group for three objects that were previously listed as exoplanet candidates: Gliese 876 b, 
HD 136118 b and HD 33636 b. Surprisingly each object has been found to belong to a different 



class: a giant planet, a brown dwarf and a M dwarf star, respectively, Benedict et al. 



(2002b), this paper, and Bean et al. (2007). These results demonstrate the importance of 



the application of complementary techniques in observing extrasolar planetary systems. 
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Table 1. Properties of HD 136118. 



ID 


HD 136118 


unit 


ref 








a 


Ljccyzviju j 


01 ■ qc; ■ Kon/f; Q7\ 
— Ul . oo . oz.oyu ^O.o/ ) 


d:m:s 


a 


/'a 


i o/i i fn ca 


mas yr ^ 


h 


/ ' 5 


ZO.O (^U. / ) 


mas yr ^ 


h 


Kabs 


19. l (0.8) 


mas 


h 


1 


Q ft Id 1 \ 
— O.O (U.l ) 




a 


Spc type 


r y V 




h 


Age 


4.8 (_ 19 ) 


Uyr 


e 


[Fe/H] 


-0.010(0.053) 


dex 


g 


[C/H] 


0.049 (0.081) 


dex 


g 


[O/H] 


0.112(0.045) 


dex 


g 


[Si/H] 


-0.042 (0.058) 


dex 


g 


[Ca/H] 


-0.057(0.062) 


dex 


g 


d 


52.3(0.6) 


pc 


h 


i> sin i 


7.33(0.5) 


ms — 1 


f 


Prot 


12.2 


day 


c 


T eff 


6097 (44) 


K 


f 




4.16(0.09) 


cm s~ 2 


b 


Af» 


1.24 (0.07) 


M Q 


c 


ft, 


1.58 (0.11) 


R 


b 


BC 


0.01 (0.03) 


mag 


b 


Afv 


3.34 


mag 


d 




7.432 


mag 


d 


V 


6.945 


mag 


d 




6.630 


mag 


d 


J 


5.934 


mag 


d 


ii 


5.693 


mag 


d 



^Perryman~| ( |199T| > 
t jAllende Prieto &: Lambert | l [ly99| ) 
IFischer et al. 1 1 |2002) 
^Zacharias et al. 1 12004 1 



e Age value and limits derived from isochrone method 
jSaffe et al. |2005| > 



Butler et al. 



(2006) 



^Gonzalez fc Laws | ( [2007] ) 
h This paper 
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Table 2. HET relative radial velocities for HD 136118. 



HJD - 2450000 RV (m/s) ± error 



3472.831 


432.6 


4.1 


3482.881 


421.4 


3.5 


3527.763 


407.8 


4.8 


3544.727 


392.4 


4.1 


3575.630 


394.6 


4.6 


3755.051 


319.0 


10.6 


3757.041 


320.5 


8.4 


3765.026 


312.2 


8.8 


3766.026 


313.0 


8.8 


3767.020 


321.9 


7.7 


3769.011 


321.9 


8.2 


3787.982 


329.2 


7.5 


3808.904 


319.1 


6.9 


3809.909 


322.5 


7.8 


3815.886 


342.9 


7.3 


3816.898 


334.4 


7.3 


3816.965 


338.7 


7.7 


3818.873 


324.5 


8.0 


3820.897 


337.5 


9.2 


3832.840 


331.2 


6.9 


3835.853 


333.6 


6.9 


3836.858 


334.6 


10.2 


3840.895 


321.5 


6.1 


3844.909 


328.9 


5.6 


3866.774 


335.1 


5.4 


3867.754 


332.1 


4.3 


3877.724 


328.3 


4.1 


3880.810 


339.0 


5.0 


3883.778 


329.3 


3.7 


3888.700 


330.7 


4.2 


3890.679 


333.3 


4.9 


3891.682 


333.7 


4.6 


3892.689 


329.5 


4.7 


3893.768 


325.2 


4.4 


3895.744 


341.4 


4.4 


3897.749 


332.9 


4.6 


3898.678 


339.1 


4.5 


3901.740 


336.0 


4.3 


3905.734 


341.8 


5.7 


3911.730 


333.9 


5.8 


3917.689 


345.5 


4.8 


3938.639 


341.2 


16.4 


3937.648 


350.4 


4.6 


3939.631 


338.5 


4.7 


4129.036 


598.0 


9.3 


4131.023 


587.5 


9.0 
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Table 2 — Continued 



HJD - 2450000 


RV (m/s) 


± error 


4135.035 


601.4 


10.3 


A A A A III > L i 

4144.998 


619.0 


8.0 


A A r- A (11(1 

4164.019 


660.3 


7.7 


4176.992 


696.4 


7.3 


4180.889 


698.6 


6.1 


4186.887 


711.6 


5.5 


4190.869 


711.2 


6.5 


4191.864 


713.7 


5.9 


4211.816 


734.9 


5.8 


4221.789 


744.9 


5.7 


4253.699 


745.3 


4.0 


4282.631 


735.8 


4.9 


4556.884 


513.6 


9.6 


4565.914 


510.6 


8.7 


4574.895 


509.8 


7.2 


4580.893 


491.0 


7.2 


4606.803 


474.3 


6.7 



Table 3. Log of astrometric observations. 



Epoch 


Date 


N ohs 


EST Roll 


1 


2005/Jun/15 


4 


58.00 


2 


2005/Jun/16 


4 


58.00 


3 


2005/Jun/17 


4 


58.00 


4 


2005/Jun/18 


4 


58.00 


5 


2005/Jun/19 


4 


58.00 


6 


2005/Jun/24 


4 


59.10 


7 


2006/Mar/02 


4 


261.00 


8 


2006/Mar/10 


4 


264.17 


9 


2006/Mar/15 


4 


266.00 


10 


2006/Mar/22 


4 


269.15 


11 


2006/Apr/03 


4 


274.00 


12 


2006/Apr/07 


4 


280.41 


13 


2007/Mar/03 


4 


261.00 


14 


2007/Mar/09 


4 


263.67 


15 


2007/Mar/15 


4 


266.00 


16 


2007/Mar/24 


4 


269.74 


17 


2007/Apr/01 


4 


274.00 


18 


2007/Apr/08 


4 


280.70 
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Table 4. Classification spectra and photometric information for HD136118 and the 

astrometric reference stars. 



Star 


Sp Ty 


V 


B - V 


My 


Av 


HD 136118 


F9V 


6.93 


0.55 


3.34 


0.0 


REF-14 


KOV 


13.95 


0.86 


5.88 


0.12 


REF-16 


GOV 


12.46 


0.73 


4.2 


0.45 


REF-17 


K0.5III 


13.55 


1.13 


0.65 


0.21 
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Table 6. HD 136118: RV orbital parameters for a 1-companion model. 



Parameter 


T T T ~\ -\ O f -1 "I i~) 1 

HD 136118b 


P [days] 


1188.0 ±2.0 


T[JD] 


2450614.7 ±6.3 


e 


0.34 ±0.01 


"[°] 


317.1 ± 1.3 


K [m/s] 


210.9 ± 1.6 


Msin? [Mj] 


12.00 ±0.47 


r^ ST [m/ s] 


485.7 ± 1.8 


^Licki [m/s] 


-1.1 ±3.8 


rLi cfc2 [m/s] 


-15.3 ± 11.9 



Table 7. Statistical Quantities (SQ) for RV residuals from 1-companion orbit fit model. 



SQ 


LICK 


HET 


ALL 


Mean 


-1.21 


-0.04 


-0.17 


Median 


0.31 


0.73 


0.27 


a [m/s] 


15.89 


8.41 


12.03 


x 2 


21.29 


92.42 


117.68 


DOF 


24 


57 


88 




0.89 


1.62 


1.34 



Table 8. Statistical Quantities (SQ) for RV residuals from 2-companion orbit fit model. 



SQ 


LICK 


HET 


ALL 


Mean 


-0.61 


0.62 


0.34 


Median 


-1.38 


0.31 


0.15 


a [m/s] 


15.4 


5.6 


10.3 


x 2 


19.9 


54.2 


76.6 


DOF 


19 


52 


83 




1.05 


1.04 


0.92 
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Table 9. HD 136118: RV Orbital Parameters for a 2-companion model. 



Parameter 


HD 136118 b 


Nuisance Orbit 


P [days] 


1191 ± 1.7 


255.3 ± 1.6 


T[JD] 


2450611 ±5 


2453761 ± 7 


K [m/s] 


215.8 ± 1.3 


11.2 ±1.9 


e 


0.353 ±0.008 


0.50 ±0.11 




316.4 ±0.9 


198 ± 16 


Msini [Mj] 


12.23 ±0.47 


0.35 ±0.07 


a sin i [AU] 


2.35 ±0.16 


0.85 ±0.25 


T HET [m/s] 


485. 


5± 1.1 


r L i C fci [m/ s] 


-1.4 ±3.3 


F L ick2 [m/ s] 


-11.4 ±9.7 



Table 10. HD 136118 b parameters obtained from simultaneous RV and astrometry 

solution. 



Parameter HD 136118 b 

P[days] 1190.9 ±1.2 

T[JD] 2450610.5 ±3.7 

e 0.352 ± 0.006 

lu[°] 316.4 ±0.6 

K [m/s] 215.99 ± 0.92 

a, [mas] 1.45 ±0.25 
163.1 ±3.0 

Q [°] 285 ± 10 

a b [AU] 2.36 ± 0.05 

M b a [Mj] 421}J 

M b a [M ] 0.0411°;™ 



a It is assumed M s = 1.24±0.07 M, 
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Fig. 1. — FGS-lr photometry of HD136118. Magnitude variation is relative to the mean 
magnitude, V=6.93. Dashed lines show the amplitude of variation possible from a (single) 
spot filling factor of 6%, the spot filling factor required to produce the observed RV variation 
from HD136118 b (ISaar & Donahue 111997b. 
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Fig. 2.— HD 136118 RV HET (filled circles) and Lick (open circles) data and the best fit 
model (solid line). 
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Fig. 3. — Histogram of residuals from 1-companion model for 3 datasets: all combined (ALL) 
(top panel), Lick (middle panel) and HET (bottom panel). 
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Fig. 4. — Lomb-Scargle periodogram for RV residuals from 1-companion model and using 
all datasets combined (solid line). The thresholds for false alarm probability of 1% and 10% 
are plotted in dotted lines. Above it is shown the AIPD (see text) for the same dataset 
(dashed line). The AIPD here is multiplied by a factor 20 and shifted for the sake of better 
visualization. 
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Fig. 5. — Lomb-Scargle periodogram for HET RV residuals from 1-companion model (solid 
line). The thresholds for false alarm probability of 1% and 10% are plotted in dotted lines. 
Above it is shown the AIPD (see text) for the same dataset (dashed line). The AIPD here 
is multiplied by a factor 20 and shifted for the sake of better visualization. 
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Fig. 6. — x 2 map from a 2-companion fit model for HET and Lick data. The grid resolution 
is about 0.3 day (1100 x 90 points). Contour lines show four different levels of x 2 ■ The best 
fit solution has the lowest value at x 2 — 76.6. 
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Fig. 7. — Phase diagram folded with period 1190.8 days. RV HET (filled circles) and Lick 
(open circles) data subtracted the "nuisance orbit" model. Solid line shows the best fit 
HD 136118 b orbit model. Residuals from the 2-companion model is plotted in the bottom 
panel. 
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Fig. 8. — RV HET (filled circles) and Lick (open circles) residuals from HD 131168 b orbit 
model and the best fit "nuisance orbit" model plotted in the phase diagram folded with 
period 255.32 days. Residuals from the 2-companion model is plotted in the bottom panel. 
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Fig. 9. — Open circles are the astrometric residuals of HD 136118 in Ax (upper panel) and 
Ay (bottom panel) versus time. Filled circles with error bars represent 3 epochs. These were 
obtained through the median and standard deviation of each grouped data. The perturbation 
orbit fit model is plotted in dashed lines. 
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Fig. 10. — Filled circles with error bars are the median and standard deviation of three 
groups of astrometric residuals of HD 136118, representing three different epochs. Dashed 
line is the fit model of the apparent perturbation orbit of HD 136118. Open circles are the 
positions calculated from the fit model, each of which is connected by a solid line to its 
respective observed epoch. The open square shows the predicted position of the periastron 
passage. 
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Fig. 11. — Histogram of Ax (left panel) and Ay (right panel) astrometric residuals for all 
reference stars. Dashed lines show the fit model for a normal distribution. 
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Fig. 12. — Predicted emission spectrum of HD 136118 (dash-dotted line), emission/reflection 
spectrum of HD 136118 b and the flux ratio between the brown dwarf and the parent star 
for a 5 Gyr and 1 Gyr system age, as indicated in the legend. 



